Stretched exponential relaxation in three dimensional short-range spin glass 
Cuo.5Coo.5Cl2-FeCl3 graphite bi-intercalation compound 



Itsuko S. Suzukj^] and Masatsugu Suzukjt] 
Department of Physics, State University of New York at Binghamton, Binghamton, New York 13902-6000 

(Dated: October 24, 2008) 

Cuo.5Coo.5Cl2-FeCl3 graphite bi-intercalation compound is a three-dimensional short-range spin 
glass with a spin freezing temperature Tsg (= 3.92 ± 0.11 K). The time evolution of the zero-field 
cooled magnetization MzFc{t) has been measured under various combinations of wait time (tm), 
temperature (T), temperature-shift (AT), and magnetic field {H). The relaxation rate SzFcit) 
[= (1/ H)dMzFc{t)/A\nt\ shows a peak at a peak time tcr- The shape of SzFc{t) in the vicinity of 
tcr is well described by stretched exponential relaxation (SER). The SER exponent b and the SER 
relaxation time tser are determined as a function of t^, T, H, and AT. The value of 6 at T = Tsg 
is nearly equal to 0.3. There is a correlation between tser and 1/6, irrespective of the values of 
tw, T, H, and AT. These features can be well explained in terms of a simple relaxation model for 
glassy dynamics. 

PACS numbers: 75.50.Lk, 75.40.Gb, 75.30.Kz 



I. INTRODUCTION 

Recently the aging phenomena have been the subject 
of many experimental studies on slow d ynamics in a va- 
riety of spin glass (SG) systems .^^^^ISEl Typically it can 
be observed in the evolution of a zero-field cooled (ZFC) 
magnetization MzFc{t) with time t after the ZFC ag- 
ing protocol for a wait time t^,. The aging behavior 
can b e und erstood based on a phenomelogical domain 
model.'^I^In this picture, the aging involves the growth 
of the domain (denoted by R) during the ZFC protocol 
for a wait time tw The domain grows with time. The size 
of the domain R becomes equal to R{tw) after the wait 
time tu,. Through this process, only the relaxation time 
T, which is nearly equal to t^j, can be selected. At t = 
just after the ZFC protocol, a magnetic field is turned on. 
Then the ZFC magnetization is measured as a function 
of the observation time t. The size of the domain (i?) 
remains constant R(t^) for < t < tyj. In contrast, the 
probing length scale {L) of the domain grows with the 
time t, starting from t = in a similar way such that the 
domain (size R) grows for the wait time during the 
ZFC protocol. The equilibrium dynamics is probed since 
L < R{tw) for < t < tnj, while the non-equilibrium 
dynamics is probed for t > tw 

An usual way to describe the slow relaxation of the 
ZFC magnetization is to postulate a statistical distri- 
bution of the relaxation times and to assume additive 
contributions. According to Lundgren et alP^ the ZFC 
magnetization MzFc{twit) is described by a sum of ex- 
ponential decay exp(— t/r) with the relaxation time r 
multiplied by the density of relaxation times g{tw, r). 



H 



[MzFcit^.t) - Mo] ^ q{t^,t) 



= g{t^,T)exp{--)dT, (1) 

where H is the magnitude of an external magnetic field. 



Mq is the ZFC magnetization at t — 0, and tq is a micro- 
scopic relaxation time (tq ~ 10"^^ sec). The relaxation 
rate SzFc(tw,t) can be defined as 



SzFc{tw, t) 



1 dMzFc{iw,t) _ dq{tu,,t) 
H dint ~ dint 

g{t^,T)- exp{-~-)dT. 

T T 



(2) 



Here it is noted that a part of the integrand expressed 
by fi^) — (1/a;) exp(— l/x) has a maximum at x = 1.0, 
where x — r/t. Using an assumption that /(r/tjis ap- 



proximated by a Dirac-delta function S(t~ t)]} 
get 



mm 



we 



g{tw,T)d{t - T)dT = g{tw,t). (3) 



Experimentally it is well known that SzFcititw) has 
a relatively flat peak centered around t — t^- This 
implies that the density of the relaxation time g(tw,T) 
also exhibits a broad peak around t = because of 
g{tw,T) = SzFc{tw:T), reflecting the glassy state. Here 
the label tu, is used for the notation of MzFc(tunt) and 
g{tu},t), in order to emphasize that each relaxation rate 
SzFc{tw,t) [~ g{tw,t)] after the ZFC protocol for the 
wait time represents the dominant feature of the ag- 
ing dynamics for a specific domain which grows for the 
wait time during the ZFC protocol. A set of data on 
SzFc{tw,t) as a function of t, for various (=10^ — 10^ 
sec) provides an information on the aging behavior of do- 
mains whose size depends on the wait time. For short t^, 
one can get the aging behavior for small domains, while 
for long tuj, one can get the aging behavior for large do- 
mains. The relaxation mechanism for the small size do- 
mains is considered to be rather different from that for 
the large size domain. There is a crossover between the 
thermal-equilibrium dynamics inside the domains and 
the non-equilibrium dynamics in domain walls. 
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FIG. 1: (Color online)Plot of SzFc{t)/Smax vs ^ = t/tc- 
which is derived from the combination of Eqs.([5]l, ([6|, ([8|, 
and Q. a = 0.035. b is changed as a parameter; b = 0.20, 
0.25, 0.30, 0.35, and 0.40. tcr is the peak time of SzFc{t). 
b>Aa^ 0.14. 



Cuo.5Coo.5Cl2-FeCl3 graphite bi-intercalation com- 
pound (GBIC) magnetically behaves like a 3D short- 
ranged SG. This compound undergoes a SG transition at 
TsG = 3.92 ifc0.11 K in the absence of H. In our previous 
papersP^we have undertaken an extensive study on the 
aging behavior of the SG phase of our system from the 
time dependence of MzFcit) under various kinds of con- 
ditions where tw, T, H and AT (the T-shift) are changed 
as parameters (see Sec.|lV]for the detail of experimental 
procedure). The relaxation rate SzFc(t) shows a peak 
at a peak time tcr- The value of t^r depends on the pa- 
rameters tw, T, H, and AT. 



In the present paper, we show that the t dependence of 
SzFc{t) is well described by the SER form in the vicin- 
ity of i ~ tcr, irrespective of t^j, T, H, and AT. The 
peak time tcr is equal to the SER relaxation time tser- 
The least squares fit of the data in the vicinity of i = tcr 
to the SER form yields the SER relaxation time tser 
and the SER exponent b, and the SER maximum S^^. 
Our results are summarized by the following two fea- 
tures, (i) the SER exponent b increases with increasing 
T. (ii) There is a strong correlation between tser and 
1 /b. These features are seen in many systems other than 
SG's, and is considered a signature of the glassy relax- 
ation. We show that these features can be well expla ined 
in terms of a simple model of glassy relaxation .112111] 



II. RELAXATION RATE S'zFc(i,.,i) 
A. General form 

Theoretically and experimentally it has been accepted 
that the the time variation of MzFcitwTt) may be de- 
scribed bjff^ 



H 



[MzFc{t^o,t) - Mo] = -At-"exp[-(i/T)'], (4) 



around t tw, where a (a > 0) is called the pre- factor 
exponent and 6 (0 < 6 < 1) is called the SER exponent, A 
is a constant, and t^ is the wait time. Then the relaxation 
rate SzFcitio,t) can be derived as 



1 dMzFc{tw,t) 
SzFc{t^,t)-- — 

= At-''cxp[-{t/T)''][a + b{t/T)% 



SzFc{tw,t) has a local maximum [5„ 
peak time tcr is given by 



(5) 

] aX t — tcr- The 



tr 



1 a ^b3{b-Aa) 
4 5 262 ) ' 



(6) 



under the condition that 6 > 4a. When a = 0, we have 



tcr 
T 



1. 



(7) 



The maximum Smax is also given by 



= _2-i+«/''[62 + ^b^b-4a)] 



-2ab + b^ + ^63(6 - 4a) 



/b 



,-2ab+b^ + ^b^b-4a). 
xexp[ ^ ]r 



(8) 



Figure [T] shows the normalized relaxation rate defined 
by SzFC / Smax as a function of the normalized time f 
defined by 



(9) 



where a = 0.035 (which is appropriate for our system, 
see Sec. IV A| and b is changed as a parameter. The 
normalized relaxation rate Sz fg I Smax has a peak at ^ = 
1 . The width of the peak in Sz fc I Smax vs ^ becomes 
narrow when b becomes increases. 



B. SER form for the least-squares fitting 

In the case of intermediate t^, one can find the time 
dependence of the relaxation rate for the relaxation of the 
sufficiently large domain size. Since the power form t"" 
is a slowly varying function of t because of small value 
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a (= 0.035) in our system, MzFc{tw,t) is approximated 
by the pure SER form, 

MzFc{t^,t) = Mo - AeM-i — )% (10) 

TSER 

in the narrow time regime near 1 = 1^, where the pref- 
actor a is equal to 0, 6 is the SER exponent and tser is 
the SER relaxation time and is nearly close to i^,. Then 
the relaxation rate SzFc{tw,t) is expressed by 

SzFc{t^,t) = eS^f^{ — )''eM-i—)% (H) 
Tser tser 

where e is the base of natural logarithmic. This func- 
tion is almost symmetric with respect to the axis t/rsER 
= 1 in the logarithmic scale of the t/tsER. axis. Since 
the density of the relaxation time g(tw, t) is nearly equal 
to SzFci^w^T), g{tw,T) has a broad peak in the very 
vicinity of r = t„,. In other words, a set of data on 
SzFc{tw,t) as a function of t, for various provides 
an information on the aging behavior of only domains 
whose size depends on the wait time. Then the expo- 
nent b should be determined from the least-squares fit of 
the data of Sz fc {tw ,t) vs < in the very vicinity of t-^j . 
To this end, in the present analysis, typically we use the 
data for the limited time regime 0.1 < < 10. Note 
that technically it very difficult to determine both a and 
b from the least-squares fit of the data of SzFcitw, t) vs 
t in the very vicinity of ■ One of the reason is that b is 
strongly dependent on the slight change of the exponent 
a. 



III. EXPERIMENTAL PROCEDURE 

The detail of sample characterization and sample 
preparation of Cuo.5Coo.5Cl2-FeCl3 GBIC was provided 
in our previous papers.^ ^ ^ The DC magnetization was 
measured using a SQUID magnetometer (Quantum De- 
sign, MPMS XL-5) with an ultra low field capability op- 
tion. The remnant magnetic field was reduced to zero 
field (exactly less than 3 mOe) at 298 K. The time (t) 
dependence of the zero-field cooled (ZFC) magnetization 
(MzFc) was measured. The following ZFC aging pro- 
tocol was carried out before the measurement. First the 
sample was annealed at 50 K for 1.2 x 10^ sec in the 
absence of H. Then the system was quenched from 50 
K to T (< Tsg)- It was aged at T for a wait time 
(typically t^, = 2.0 x 10^ - 3.0 x 10^ sec). After the 
wait time, an external magnetic field H is applied along 
any direction perpendicular to the c-axis at t = 0. The 
measurements of Mzfc vs t were made under the vari- 
ous conditions. The detail of the experimental procedure 
will be described in Sec. IIVI 



IV. RESULT 
A. Determination of the pre-factor exponent a 

In the early stage of the time evolution (i ~ 0, t^, — 0, 
but still t <C tw), the ZFC magnetization can be well 
described by a power-law form: 

^[MzFc{t^,t) - Mo] = -At-\ (12) 

The corresponding relaxation rate is obtained as 

SzFcit^,t)^t-''. (13) 

The relaxation rate monotonically decreases with increas- 
ing t in the early stage. In our measurement using SQUID 
magnetometer, it is very difficult to measure the time de- 
pendence of Mz FC {tw , t) for very short ■ No rehable 
data can be taken for small since it takes relatively 
long time until for the temperature to become stable be- 
fore the ZFC process starts. There is some uncertainty 
for the definition of t^. In the early stage (t = 0), in turn, 
we use the frequency dependence of the absorption of the 
AC magnetic susceptibility to determine the value of a. 
The absorption x"(w) of the AC magnetic susceptibility 
is related to the relaxation rate through a so-called tt/2 
law,E2l 

X"(^) = -'^SzFc{t^u - l/^: t = 0) « (14) 

In our previous papei'^ we experimentally determine the 
value of a as 

a = 0.035 ±0.001. (15) 

In fact, it is experimentally confirmed that this exponent 
a is related to the critical exponents b, and ziy by 

a=-, (16) 
zv 

where /3 is the critical exponent of the order parameter 
(/3 = 0.36 ±0.03), z is the dynamic critical exponent, and 
V is the critical exponent of the inverse correlation length 
in the present system; zv = 10.3 ± 0.7. 

B. Aging behavior at various T 

We have measured the t dependence of MzFc{t) at 
various T after the ZFC aging protocol where H — 1 
Oe and i^, = 2.0 x 10^ sec."^^ Figure ^a.) shows the t 
dependence of SzFc{t) at various T, where H — 1 Oe 
and tuj — 2.0 x 10^ sec. The relaxation rate SzFc{t) 
shows a peak (the peak height Smax) at the peak time 
t — tcr, which shifts to the short-t side with increasing T. 
The peak time t^r drastically decreases with increasing 
T near T = Tsg- The existence of the peak at t = tcr 
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FIG. 2: (Color online) (a) t dependence of SzFc{t) at various 
T for Cuo.5Coo.5Cl2-FeCl3 GBIC. 3.0 < T < 4.7 K t„ = 
2.0 X 10^ sec. H = 1 Oe. The solid lines denote the least- 



squares fits to a SER given by Eq.(ff I for Szfc vs t in the 
vicinity of the peak time tcr- The fitting parameters S'^ax, b, 
and TSER in Eq.(|ll[) thus determined are shown in Figs.jsja), 

(b) , and (c) as a function of T. tser is the relaxation time 
for the SER. b is the SER exponent, (b) Scaling plot of the 
ratio SzFc(f) / Smax as a function of t/tcr for 3.3 < T < 4.5 
K. The values of tcr and Smax are shown in Figs. |3|a) and 

(c) , respectively. H = 1 Oe. = 2.0 x 10^ sec. 



in SzFc{t) vs t indicates that the SER plays a signifi- 
cant role around t = tcr- We find that SzFc{t) is well 
described by the SER form given by Eq.pT]). The least- 
squares fit of these data of SzFcit) vs t to Eck|iTJ yields 
the parameters 6, tseRi and S'^f^. Figure |2|b) shows 
the plot of SzFc{t) I Smax as a function of t/tcr at various 
T, where Smax and tcr are different for different T. 



Figures [sjja) 
Tser, b 
fits, where Tsg 



and S^jf^ , 



(b), and (c) show the T dependence of 
determined from the least-squares 
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FIG. 3: (Color online) (a) tcr vs T and tser vs T. (b) b vs 
T. (c) 5"^^^ vs T and S'^ao; vsT. H ^ 1 Oe. t„ = 2.0 x 10^ 
sec. tcr is a characteristic time at which SzFc{t) exhibits a 
peak at the fixed T. Smax is the peak height of SzFc{t) at 
t = tcr- The arrows indicates the location of Tsa (= 3.92 K). 



1 Oe. The values of tcr and Smax are also plotted as 
a function of T. Note that the peak time tcr and the 
height Smax are derived directly from the t dependence 
of SzFc{t)- We find that tser drastically decreases with 
increasing T in the vicinity of Tsg- The value of tser. at 
each T is almost the same as that of tcr at the same T. It 
should be noted that tcr (or tser) is equal to t^ &t T = 
Tsg- In Fig- (lit)) we show the T dependence of b. The 
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FIG. 4: (Color online) Plot of 1/tser and 1/tcr as a function 
of l/T. H = 1 Oe. t„ = 2.0 x 10^ sec. The solid lines denote 




the best fit of the SER form given by Eq.( 11 1 to experimental 
data for T < Tsa and T > Tsg, respectively. The arrow 
indicates the location of Tsg- 



exponent b is equal to 0.15 aX T — 3.3 K and increases 
with increasing T. The exponent h is nearly equal to 0.3 
at T = Tsg- This value of h is in good agreement with 
the value at Tsg which is predicted by Ogielski^^ from 
the Monte Carlo simulation. Similar behavior has been 
observed by Bontemps and Orbach^^ and Bontemp^^^ in 
the curve b vs T for the insulating spin glass Euo.4Sro gS 
{Tsg = 1-5 K): the exponent b is nearly equal to 0.17 at 
1.3 K. It increases with increasing T. It is equal to 0.3 
sXT — Tsg and increases to 0.4 just above Tsg- 

In Fig. ^c) we show the values of S^ax and S^^^ as 
a function of T. The value of Smax is almost equal to 
that of S'mai' at the same T. We find that Smax shows a 
broad peak at T = 3.6 K just below Tsg- The deviation 
of the data of SzFc{t) t from the SER occurs at both 
t ^ tcr and t <^ tcr- For convenience, in Fig. [2{b) we 
define characteristic times i„ (> tcr) and t; (< icrTT^tiGre 
SzFc{t) reaches a O.SSmax', ^"^(tu/tcr) — 0.6012/6 and 
lii(ti/tcr) = —0.7515/5. The times i^, and t; approach 
the peak time tcr as b is increased. It follows that the 
width of the curve SzFc{t) / Smax vs tjtcr, defined by 
ln(t„/t;) — 1.35/6, becomes narrower as b is increased 
from 6o = 4a = 0.14 to 6 = 1. The width of the curve 
{SzFc{t) / Smax VS t) becomes narrower with increasing 
T. This implies that b increases with increasing T (see 
Fig. |3]^b) for comparison). 

In Fig. |4] we show the plot of 1/tser (or 1/tcr) as a 
function of l/T. There is a drastic change in 1/tser vs 
l/T around l/T ^ 1/Tsg ^ 0.255 K"^ It is expected 
that the t dependence of tser (or tcr) is given by an 
Arrhenius law. 



with different and for T > Tsg and T < Tsg- 
The least-squares fit of our data of of 1/tser (or 1/icr ) 
yields the parameters Cq — 19. 7± 11.5 sec~-^ and — 
10.82 ± 0.56 for T > Tsg, and = 1.01 ± 0.48 sec"i 
and cl = 7.81 ± 0.46 for T < Tsg- The temperature 
corresponding to the characteristic energy barrier of the 
relaxation process is Eb = c^Tsg = 42.4 K for T > Tsg 
and 30.6 K for T < Tsg- The same form of the 1/tser 
vs l/T has been used by Hoogerbeets et aL^^^ for their 
analysis of TRM relaxation measurements of canonical 
SG systems: Ag:Mn (2.6 at. %), Ag: Mn (4.1 at. %), 
Ag:[Mn (2.6 at. %) + Sb (0.46 at. %)], and Cu:Mn (4.0 at 
%). Our value of ct is much larger than those derived by 
Hoogerbeets et al]^ {cl — 2.5). Note that in their work 
the stretched exponential was taken as representative of 
the short time {t < tw) relaxation. 



C. Aging behavior at various 

We have measured the t dependence of MzFc{t) at T 
= 3.75 K and = 5 Oe just below Tsg after the ZFC 
aging protocol, where the system was aged at T = 3.75 
K for the wait time tw This wait time t^ is varied as 
a parameter: 1.0 x 10^ < tyj < 3.0 x lO'^ sec. A least- 
squares fit of these data in the vicinity oi t — tcr to 
Eq.(ll) for t/tcr < 10 yields the parameters tser and 
b for each tyj. Figures ^a) and (b) show the plot of 
TsERi ten and 6 as a function of t^. We find that the 
exponent b decreases from 0.4 to 0.2 with increasing t^j 
from 1.0 X 10^ sec to 3 x 10'' sec. The data of b vs tyj for 
1.0 X 10^ <t^ < 3.0 X 10* sec is described by 



b = b*o-bllii{t^/t"J, 



(18) 



102 







0.40 ±0.02, and 
the value of b can 
2.0 X 10^ sec, which is 



where is chosen as 1.0 x 
bl = 0.034 ± 0.004. Using Eq.([l8l 
be estimated as 6 = 0. 
nealy equal to 4a. 

In Fig. ^c) we show the plot of 1/tser as a function 
of It is suggested by Chamberlin.E^ that tser. is 
described by a form 



1/tser = w* exp(-t^/i* ), 



(19) 



where lo* and t*, are constant to be determined. A least- 
squares fit of the data of 1/tser vs t^^ in t he limited 
i^-region 750 < i^, < 1.5 x 10"* sec to Eq.(19) yields 
(3.86 ± 0.25) X lO-'' 



750 < ! 
the parameters 



1/tser = Cq exp{-clTsG/T), 



(17) 



X lU ■* sec and 
tl = (1.13 ± 0.18) X 10* sec. The values of 1/tser for 
tw < 750 sec considerably deviates from the exponential 
tyj dependence, partly because of the initial stage of the 
aging process depending on the initial condition. 



D. Aging behavior under the T-shift 

We have measured the t dependence of MzFc{t) under 
the T-shift from the initial temperature Ti = T — AT 
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FIG. 5: (Color online) (a) t^r vs and tser vs tyj. (b) the 
SER exponent 6 vs T = 3.75 K. = 5 Oe. 1.0 x 10^ < 
< 3.0 X 10* sec. The solid line is the best fit of Eq.([l8l 
to the data of 6 vs These results are obtained from the 
measurement of MzFc{t) as a function of t after the ZFC 
cooling protocol and isothemal aging at T (= 3.75 K) for a 
wait time tw (c) 1/tser vs and 1/tcr vs t^. T — 3.75 K. 
H = 5 Oe. The solid line is the best fit of Eq.( 19 1 to the data 
of 1/tser vs t„. 
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FIG. 6: (Color onhne) (a) tcr vs AT and tser vs AT. The 
solid line denotes the best fit of Eq.( |20[ ) to the data of tser 
vs at for < AT < 0.25 K. (b) b vsAT. H = 5 Oe. U = 
3.0 X 10* sec. The results are obtained from the measurement 
of XzFc{t) as a function of t, after the ZFC cooling protocol, 
isothermal aging a,t T — Ti — Tj ~ AT. Immediately after the 
temperature is shifted to T — Tf (so-called AT shift) and the 
magnetic field is turned on, Mzfc is measured as a function 
oit. 



to the final temperature Tf = T, where Ti = 3, 3.2, 
3.4, 3.5, 3.6, and 3.9 K. In the T-shift experiment, the 
system is cooled in zero field from 50 K to a temperature 
T, = T- AT below TsG- After a wait time (= 3.0x10'' 
sec) at this temperature, immediately prior to the field 
application, the temperature is raised to Tf = 3.75 K. 
The t dependence of MzFc{t) was measured a,t H = 
5 Oe. We have already reported the t dependence of 
SzFc{t) under the T-shift in our previous paper. ^ We 
find that Szpcit) shows a peak at tcr- The peak shifts 
to the long-i side as AT is decreased. A least-squares 
fit of the data SzFc(t) vs t in the vicinity of t = tcr to 
Eq.(ll) yields the parameters tser and b. In Figs. [6]^ a) 
and (b), we show the parameters tser (also tcr) and b 
as a function of AT. The relaxation time tser decrease 
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FIG. 7: (Color online) Scaling plot of the ratio SzFc{i) / Smax as a function of t/tcr- tw = 1-0 x 10* sec. H is changed as a 
parameter, (a) T= 3.1 K, (b), 3.3 K, (c) 3.5 K, and (d) 3.65 K. 



with increasing AT, while b increases with increasing AT. 

In our previous paper,* we have shown an expression 
for tcr under the AT-shift aging process. This expres- 
sion is derived from the Monte Carlo simulation based 
on the droplet model (Takayama and Fukushima^*) un- 
der the condition that the domain size is comparable to 
the overlap length Xat- In the limit of AT --^ 0, Ihtser 
is linearly dependent on AT, 



TSER. = T-5;exp(-a^AT), 



(20) 



where a slope and a relaxation time are to be 
determined. The slope increases with increasing tw 
The curve {Iutser vs AT) for < AT < 0.25 K is 
linearly dependent on AT (see Fig. |6|a)). The fitting 
parameters are given by slope — 8.7 ± 0.5 K^^ and 
= (1.81 ±0.05) X 10" sec. 



E. Aging behavior under the //-shift 

We have measured the t dependence of MzFc{t) after 
the ZFC aging protocol which consists of (i) annealing at 



50 K for 1.2 X 10^ sec at F = 0, (ii) cooling from 50 K 
to T (< Tsg), (iii) isothemal aging at T for a wait time 
{— 1.0 X 10^ sec), and (iv) switching a field from to 
H. The time t = is a time when H is turned on. Our 
results on the dependence of SzFc{t) at the fixed T and 
H of the {H,T) plane, have been reported in our previ- 
ous paper.- The relaxation rate SzFc{t) exhibits a peak 
[the peak height Smax{T, H)] at a peak time tcr{T,H), 
which drastically shifts to the short-t side with increas- 
ing H. Figure [t] shows the scaling plot of SzFc{t)/ Smax 
as a function of t/tcr a± T — 3.1, 3.3, 3.5, and 3.65 K. 
The magnetic field is changed as a parameter. The least- 
squares fit of the data of SzFc(t) vs t in the vicinity of 
tcr {T, H) to the SER given by Eq. ( [TT) , yields the parame- 
ters TsER and b. In Figs.[8|a) we show the H dependence 
of TsER at various T below Tsg- Note that the values of 
'TSER, and tcr are almost the same at the same T and H. 
We find that tser. drastically decreases with increasing 
H at each T below Tsg- In Fig- |8jb) we show the H 
dependence of b at various T. The exponent b increases 
with increasing H at each T and reaches a value between 
0.4 and 0.5. The field H at which b is equal to 0.3, in- 
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FIG. 8: (Color online) (a) r vs Ji" at various T. i„ = 1.0 x 10* 
sec. The lines denote the best fits of Eq.(21l to the experi- 
mental data in the low H limit, (b) 6 vs H. The lines are 
guides to the eyes, = 1.0 x 10* sec. The results are ob- 
tained from the measurement of MzFc{t) as a function of 
t in the presence of H, after the ZFC cooling protocol and 
isothermal aging at T for a wait time 'At H = 0. 



creases with decreasing T from 3.75 to 2.9 K. Note that h 
decreases with increasing the coohng field He just below 
TsG for the TRM decay experiment on Ag:Mn(2.6 at %) 
-I- Sb(0.46 at %) (Hoogerbeets et aXW^ and Cu: Mn (6.0 
at %) (Chu et al.ESl). This result is different from our 
result from the ZFC magnetization relaxation measure- 
ment that h increases with increasing the applied field 
H. 

In our previous paper P we have shown an expression 
for tcr under the iJ-shift aging process. This expres- 
sion is derived from the Monte Carlo simulation based 
on the droplet model (Takayama and FukushimaP^l) un_ 
der the condition that the domain size is comparable to 
the crossover length L^- In the limit of if ^ 0. Inf^r is 
linearly dependent on H: Inter ~ —ctnH, where is 



FIG. 9: (Color online) T dependence of fitting parameters H* 
and r|f determined from the least-squares fit of the data of 
Fig.[^a) to Eq.([2l]). (a) H* vs T. (b) vs T. 



constant. We assume that tser is described by the form 



TSER = exp{~H/H*), 



(21) 



8]^ a) we 
ow TsG 



in the low-iJ limit, where H* = In Fig. 

show the H dependence of tser at various T bel 
where tw = 1.0 x 10^ sec. We find that In tser is propor- 
tional to H at low H. The slope aj^ gradually decreases 
as T is lowered for T < TgQ. The least-squares fit of 



the data of tser vs H at low H to Eq.(21) yields pa- 



rameters r|f and H* at each T. In Fig. |9J we show the 
T dependence of H* and r|j . The characteristic field H* 
decreases with increasing T below TgQ. The character- 
istic relaxation time t'^ decreases with increasing T and 
reduces to zero around T = Tsg- 



V. DISCUSSION 

We have shown that SzFc{t) exhibits a peak at a peak 
time tcr- The relaxation rate SzFc{t) is well described 
by a SER in the vicinity of . The relaxation time ts er 
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FIG. 10: (Color online) Relation ofrsER vs 1/b. All the data 
of TSER vs b obtained in the present work are plotted. The 
data fall well on a single curve given in the text, irrespective 
of the values of t^n, H, T, and AT. The data with different 
notations are obtained under different conditions (t^, T, H, 
and AT). 



obeys the Arrhenius law. The exponent b for the SER 
increases with increasing T and is nearly equal to 0.30 
at T = TsG- This value of b at Tsg is in good agree- 
ment with the prediction from the numerical simulation 
on the ±J Ising spin glasses (OgielskP^l). The exponent 
b increases with increasing T above TgQ and tends to ap- 
proach 1 (exponential type relaxation, Debye-type) well 
above Tsg- Similar conclusion has been derived by Keren 
et from the T dependence of the spin-spin dynami- 
cal autocorrelation function of the Ising SG, Feo.o5TiS2: 
b is equal to 1/3 at T = Tsg and b increases with increas- 
ing T above Tsg- The value of 6 « 1/3 at T = Tsg may 
imply that the available configuration space tends to the 
structure of a percolation fractal, so the SG transition is 
physically a percolation in phase space (Almeida et al.-^) 
Theoretically it is predicted th at th e exponent b is re- 
lated to the critical exponents h^p^^ 



6 = 



P + 1 



P + j + zv' 



(22) 



where 7 is the exponent of the susceptibility. Using our 
critical exponents {(3 = 0.36 ± 0.03. 7 = 3.5 ± 0.4 
, and zi/ = 10.3 ±0.7), the critical exponent b for 
Cuo.5Coo.5Cl2-FeCl3 GBIC can be estimated as & = 0.27, 
which is close to 0.3 at T = Tsg- Note that the crite- 
rion (6 > 4a) is satisfied for the existence of the peak of 
SzFc{t) aX t — tcr, since a = 0.035. 

Figure [To] shows the plot oirsER^s ^/b for all the data 
obtained in the present work, including Figs. [3|^a) and 
(b). Figs. [5]; a) and (b), Figs.|6|a) and (b), and Figs.[8]^a) 



FIG. 11: (Color online) Relaxation rate S{t) vs t at various 
values of ct: The parameter 1/a is proportional to T. The 
t dependence of S{t) is derived from numerical calculation 
of the differential equation given by Eq.(|24[). The relaxation 
rate S{t) is well described by the SER form defined by Eq. 1 11 1 
with b in the vicinity of the peak time. 



and (b). The data oirsER vs 1/6 for 2 < 1/6 < 3.3 (cor- 
responding to 0.3 < 6 < 0.5) fall, to within experimental 
accuracy, on a single line which is the best fit of the form. 



In rsBi? ^pI+p\ ln(l/6), 



(23) 



with pI = 0.9 ± 0.3, pI = 5.2 ± 0.5. On the other hand, 
the data of tser vs 1/6 for 3.3 < 1/6 < 7 (0.1 < 6 < 0.3) 
are rather broadly distributed near the single line. This 
result suggests that the exponent 6 is closely related to 
the SER time tser, irrespective of the val ues o f t^, 
T, H, and AT. Note that Hoogerbeets et al.ESlll have 
reported a relationship between Ihtser. and 1/6 below 
Tsg for canonical SG systems: Ag:Mn (2.6 at. %), 
Ag: Mn (4.1 at. %), Ag:[Mn (2.6 at. %) + Sb (0.46 
at. %)], and Cu:Mn (4.0 at %): They have shown that 
^T^TSER decreases with increasing 1/6. Their results are 
very different from our result. The correlation between 
TSER and 6 has been studied by Goltzer et al^ using 
Monte Carlo simulations on the Ising SG's: tser in- 
creases and 6 decreases as T is approached Tsg from the 
high-temperature side. This result is similar to our re- 
sult. 

A successful theory of the glassy relaxation should pro- 
vide a justification for two common features: (i) the SER 
form of the relaxation rate and (ii) the correlation be- 
tween Tser and 1 /6. These features are seen in many sys- 
tems other than SG's, and is considered to be a signature 
of the glassy relaxation. Here we consider a simple model 
proposed by Trachenko and Dove,'^and Trachenkd^ for 
the glassy transitions. The dynamics of the local relax- 
ation events is governed by a differential equation with a 
solution that fits well to the stretched-exponential relax- 
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FIG. 12: (Color online) Fitting parameters of the SER. (a) 6 vs 1/a, where 1/a is proportional to T. The exponent 6 is equal 
to 0.3 at 1/a = 0.175. (b) tser vs a, which denotes the Arrhenius law. (c) tser vs 1/a. The relaxation time tser increases 
with decreasing 1/a. (d) The relationship between tser and 1/6. The solid lines in (b) and (d) denote the best fitting curves 
(see the text). 



ation. The rate equation is given by 



dx{t) 
dt 



exp[— Q;x(t)] — x(t) exp(— a), 



(24) 



with an initial condition x(0) = 0, where x(t) may cor- 
respond to the ratio MzFc{t) /MzFc{t — co)- The pa- 
rameter a is proportional to Es/kBT, where Eb is the 
activation barrier, /c^ is the Boltzmann constant. The 
time t is redefined as t/to, where to is a characteristic 
time. The second term of Eq.(24l describes saturation, 
such that dx{t)/dt = as t — > oo, or x{t) 1. The 
relaxation rate S(t) is defined by 



S{t) 



dx{t) 
dint 



~ t{cxp[—ax{t)] — x{t) exp(— a)}. (25) 



Figure [TT] shows the result of numerical calculation of 
S{t), where a is changed as a parameter (a = — 10). 



The parameter (1/a) is proportional to T. The relax- 
ation rate S{t) has a peak around t — 1. This peak 
shifts to longer-i side with decreasing {1/a) (or with 
decreasing T). We find that S{t) is well described by 
the SER form given by Eq.(ll) in the vicinity of the 
peak time. The least-squares fit of the data of S{t) vs 



t to the SER form given by Eq.(lll yields the param- 
eters b and tser- Figures 12'a)shows the plot of b 
vs 1/a. The exponent b increases from b — 0.1 to 0.4 
with increasing (1/a). The exponent b is equal to 0.3 
around 1/a = 1/asG = 0.175 corresponding to Tsg 
{asG = 5.71). Figure [l2|b) shows the relaxation time 
TSER as a function of a (oc 1/T), showing the Arrhe- 
nius form. The relaxation time tser for asa < a < 10 
is well described by the form, Xyitser = d*Q + d\a with 
dl ^ -0.863 ± 0.006 and dl = 0.397 ± 0.002. The pa- 



rameter c{ for the Arrhenius law Eq.(17l is related to 



by c{ — d\asG = 2.26. This value of C2 is close to 
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that reported by Hoogerbeets et all^^ = 2.5). Fig- 



ure 12 ^c) shows the relaxation time tser as a function 
of ija ((X T). The relaxation time tser drastically in- 
creases with decreasing 1/a below 1/a = 1/asG- Figure 



12 'd) shows the relationship between tser. and 1/5. The 
relaxation time tser. is uniquely determined from the 
value of 1 /6 by the relation 



l^TsER = qo+qt ln(l/&), 



(26) 



for 3.3 < l/b < 10, where = -0.35 ± 0.02 and ql = 
1.515 ± 0.01. Our experimental value of pi (= 5.2 ± 0.5) 
is much larger than the theoretical value of ql , suggest- 
ing the incompleteness of the model. In spite of such 
difference it can be concluded from this model that er 
increases with increasing 1 /b. This is the same conclusion 
derived by Ngai and Tsang.1221 In summary, the features 
of parameters of the SER, can be explained in terms of 
the simple model of the glassy dynamics.f^^ 



time dependence of the relaxation rate SzFc{t)- The 
relaxation rate SzFc{t) is described by a SER form in 
the vicinity of t = t^r- There is a correlation between the 
exponent r and l/b of the SER, irrespective of the val- 
ues of tf^, T, H, and AT. The exponent b is equal to 0.3 
at r = TsG- The exponent b increases with increasing 
T. The relaxation time tser obeys the Arrhenius law. 
These features, which are a signature of the glass relax- 
ation of many systems, can be well explained in terms of 
the simple relaxation model. 
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